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Allometric scaling is one of the most pervasive laws in biology. Its origin, however, is still a 
matter of dispute. Recent studies have established that maximum metabolic rate scales with an 
exponent larger than that found for basal metabolism. This unpredicted result sets a challenge that 
can decide which of the concurrent hypotheses is the correct theory. Here we show that both scaling 
laws can be deduced from a single network model. Besides the 3/4- law for basal metabolism, the 
model predicts that maximum metabolic rate scales as M^^'' , maximum heart rate as M~^'^^, and 
muscular capillary density as Af in agreement with data. 

PACS numbers: 87.10.+e,87.23.-n 



Metabolic rate B and body mass M are connected 
by the scaling relation B — aAl'', where a is a con- 
stant and b is the allometric exponent. Kleiber's law 
(6 = 3/4) characterizes the basal metabolism of almost 
all organisms |^ Q, including mammals and birds. 
The predominance of quarter-power scaling in biology 
is experimentally evident |^ |^ Q and is supported by 
recent theories based on the design of resource distri- 
bution networks j5i,i^]. Maximum metabolic rate, how- 
ever, scales with an exponent somewhat larger than 3/4 
(6 « 0.86) mmiilialll. This result raises the ques- 
tion whether the scaling of maximal metabolic rate is 
governed by mechanisms different from those determin- 
ing basal metabolism 0, |M • Here we reformulate the 
ideas of West, Brown and Enquist to show that both 
scaling laws can be deduced from a single network model. 
Although it was recently demonstrated that the exponent 
3/4 cannot be consistently derived from the assumptions 
of West et al. 0,0], we show that, using a different set of 
assumptions and the appropriate impedances, the model 
correctly predicts the 3/4-law for basal metabolism. 

The basal and maximum metabolic rates set the limits 
of the energetic range of endothermic animals. Consider- 
able effort has recently been invested to understand the 
scaling of these variables "i, 's', "e", 7, Tllil [H El El 
Es. 16j. In particular, maximal performance is very infor- 
mative about animal design since during strenuous exer- 
cise the circulatory and respiratory systems are stressed 
to their uttermost. Recent studies have established that 
maximal metabolic rate (MMR) scales with an exponent 
larger than that found for basal metabolism, suggest- 
ing that the scaling of MMR is determined by mech- 
anisms different from those controlling basal metabolic 
rate (BMR) HHH. We show that this difference in the 
allometric exponent is a consequence of dynamic adap- 
tations in the resource supply network, which are needed 
to cope with the distinct demands for oxygen and nutri- 
ents. Simple dimensional analysis linking BMR to heat 
loss through the surface (claimed to support a 2/3-law for 
basal metabolism) 12] and simple geometric arguments 
are unable to explain why different metabolic states 
have different allometric exponents. 



West, Brown and Enquist (WBE) proposed that the 
ubiquity of the 3/4-law lies in the resource distribution 
network common to all organisms |^. In order to esti- 
mate the exponent b they make three assumptions: (i) 
the network branchs hierarchically to supply all parts of 
the organism; (ii) the terminal units, such as the cap- 
illares of the circulatory system, are body size indepen- 
dent; (iii) natural selection has optimized the network 
so that the energy required to distribute resources is 
minimal, or, equivalently, the total hydrodynamic re- 
sistance (impedance) of the network is minimal. From 
these assumptions they conclude that (a) the network is 
self-similar (fractal) ; (b) the network has area-preserving 
branching in the first generations; (c) blood volume Vf, 
is proportional to M; and (d) metabolic rate scales with 
6 = 3/4. Some matematical errors, however, have been 
found in the derivation of these results E3- It turns out 
that none of these conclusions can be derived from the 
above hypotheses. Nevertheless, the ideas of WBE are 
simple and elegant, and will be shown here to predict 
correctly the allometric exponents of BMR and MMR 
when a different set of assumptions is adopted and the 
correct impedances are used. 

The impedance is not minimal for pulsatile flow in an 
elastic tube, nor for an area-preserving branching in a 
rigid tube. We therefore reformulate the third hypothe- 
sis to state that (iii) the total impedance is minimal with 
respect to the radius, length, thickness and elastic modu- 
lus of the vessels. In other words, the impedance is min- 
imal with respect to all variables, except the branching 
ratio n. In addition, we assume that (iv) blood volume 
is proportional to M; (v) the network is fractal; and (vi) 
the mass of the arteries is constant for a given M or, 
alternatively, (vi) the total elastic energy stored in the 
vessels is proportional to M. Hypotheses (iv) and (v) 
are supported empirically [3,HE3iE3); assumption (vi) 
is discussed below. 

Let us consider some aspects of the cardiovascular sys- 
tem E3 and their relation to the model. For clarity, we 
use the same notation as West et al. @ . The circulatory 
network consists of elastic vessels: the aorta, large ar- 
teries, muscular arteries, arterioles and capillaries. Each 
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vessel is described by a label k related to the branching 
level, with k — corresponding to the aorta and k — N to 
the capillaries. At each level, each vessel branches into 
Uk smaller vessels, so that the total number of vessels 
Nk at level k is Nk — nQUi . . .rik-i- In particular, the 
number of capillaries is Nc — Nn = noni . . . hm-i- After 
branching, a vessel has a smaller length (Ik+i < Ik) and a 
smaller radius {rk+i < rk) and the network can be char- 
acterized by the ratios jk = h+i/lk and (3k = rk+i/rk- 
The heart is taken into account by considering an oscilla- 
tory pressure = pQexp[i{ujt~ Kz)], where po is the 
amplitude of the pulse, t is time, z is the distance along 
the vessel, K = Ki + 1X2 is the complex wavenumber 
and uj = 27r/ is the angular frequency. The velocity of 
the pressure pulse is given by c = uj/Ki = A/, with A 
the wavelength and / the frequency; K2 describes the 
damping of the pulse. 

The thickness hk and the Young's modulus Ek of the 
vessels also depend on the level k. The vessel wall con- 
sists, in general, of three layers: the intima, a monolayer 
of endothelial cells, the media, consisting of muscle layers, 
and the adventitia, which is composed of connective tis- 
sue. The elastic properties of the vessels vary principally 
due to the variable constitution of the media. Large ar- 
teries are distentensible and have a smaller E than mus- 
cular ones. Arterioles have a relatively large thickness 
and small changes in muscular contraction can induce 
radius changes, thus controlling flow to organs. The cap- 
illaries have no adventitia or media. Their function is the 
exchange between blood and cells. 

The circulatory system of endothcrmic animals is a 
dynamical network which is adjusted according to the 
metabolic state. In humans, only 20-25% of the muscular 
capillaries are perfused during rest, while all are perfused 
during exercise. This increase is due to the opening of 
terminal arterioles and to the rise in perfusion pressure at 
the capillary level. (Note that this implies that muscular 
capillary density is linked to MMR, instead of BMR.) At 
rest blood flow to skeletal muscles is only a small frac- 
tion of the total (15% in humans), but during strenuous 
activities it increases to a large fraction (over 90% in 
mammals) Q). Circulating blood volume is also larger 
during MMR, since the veins function as a reservoir dur- 
ing BMR. 

In short, the transition from resting to maximum ac- 
tivity can be described as follows: (a) the heart increases 
its rate and output; (b) the mean arterial pressure and 
peripheral extramuscular resistance increase due to the 
sympathetic vasoconstriction; (c) circulating blood vol- 
ume increases due to vasoconstriction of the veins; (d) 
extramuscular flow remains essentially constant, some- 
what reduced in some organs but increased in others; and 
(e) total flow and muscular flow increase, with all muscu- 
lar capillaries activated. The increase in heart rate and 
in circulating blood volume during maximal performance 
highlights the elastic nature of the vessels and leads to 
distinct allometric exponents for BMR and MMR. 

The exponent b of the metabolic rate for a network 



with rife = n is obtained as follows. The volume rate of 
flow in a tube is Qk = nr'^Uk, where Uk is the average 
velocity. Since the fluid is conserved we have 

Qo = NkQk = N^nrluc . (1) 

Because the quantities related to the capillaries (re, 
and Uc) are size invariant (second hypothesis), we obtain 
Qo ~ Nc- The fluid tranports all the nutrients and oxy- 
gen for the metabolism, so it is natural to suppose that 
B ^ Qo- Therefore we have that B ^ Nc ^ ^ M''. 
From hypothesis (iv), we know that M is related to the 
network quantities through the total blood volume Vb- 
Evaluating Vb explicitly in terms of 7, /3 and capillary 
quantities yields 
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and we obtain 



ln(n) 



(2) 



(3) 



In order to evaluate the exponent b we need 7 and (3. 

Let us now consider a network characterized by the 
variables {lk,rk,nk,hk, Ek}. The first hypothesis can 
be described by Nhll — Nk+il^^i, the so called space- 
filling condition which fixes the length fraction as 

— 1/3 

Ik — rif, . The ratio comes from minimization of the 
total impedance Z (third hypothesis). During maximum 
performance, the elastic nature of the vessels is important 
because circulating blood volume and volume flow take 
their largest values. Moreover, heart rate is maximum 
and the wavelength A of the pressure pulse is minimum, 
so that we must evaluate a local impedance. Analyzing 
the oscillatory motion of a viscous liquid in a thin-walled 
elastic tube, Morgan and Kiely J^J and Womcrsley [2(tI | 
showed that the local impedance Zk, defined by 



Q{z,t) 



piz,t) 



— ^p(z,t) 



(4) 



where Q{z,t) and p{z,t) are the volume flow and the 
pressure wave at position z and time t, is given by 



Zk = 



5/2 



(5) 



Here A — Pq^^ /y/2Tr with po the blood density. This ex- 
pression was obtained by solving the Navier-Stokes equa- 
tion for the liquid coupled to the Navier equation for the 
vessel wall, in the linearized long wave approximation for 
large r. Under these conditions, valid for large elastic 
arteries, the pressure wave is not attenuated {K2 = 0), 
propagating with velocity c = u)/Ki = {Ekhk/^poTkY^"^ . 

We want to minimize the power dissipated in the 
circuit (cardiac output) W = ZQ^. Since B ^ Qo, 
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TABLE L AUometric exponent b describing the dependence 
of a variable Y on body mass M (F ~ M''). 



Variable Exponent 





Predicted 


Observed 


Ref. 


Basal metabolic rate 


3/4 = 0.75 


0.737 ± 0.026 


4 


Maximum metabolic rate 


6/7 0.86 


0.828 ±0.070 


4 






0.88 ±0.02 


7 






0.872 ±0.029 


8 






0.86 


22 






0.85 


23 


Muscular capillary density 


-1/7 f» 0.14 


-0.07 to -0.21 


24 


Heart rate 


-1/7 ^ 0.14 


-0.17 ±0.02 


25 






-0.16 ±0.02 


26 






-0.15 


27 



where 77 is the viscosity. Using their results, we obtain 
the impedance 

^fe = 1/2 7/2 ' 

where C = 8/(V57r). Womersley's analysis leads 
to a similar result. We may then minimize the total 
impedance subject to the same restrictions as before, and 
within the same approximation {Ek = E) . We readily ob- 

— 1/3 

tain (3k = ■ order to determine the exponent 6, 

we employ Uk — n and a steplike /3, with [5 — n^'^l^'^ for 
k < k* and (3 = n~^^^ for k > k* , where k* is the level 
at which the two impedances are comparable. Following 
West et al. Q , we find that Vb is dominated by the large 
vessels, implying that the MMR allometric exponent is 
6 = 6/7 « 0.86. 

So far we have assumed the elastic modulus to be in- 
dependent of k. To take the variation of Ek with level 
k into account, we need a restriction that depends on 
Ek- First, we note that the circumferential tension in 
the tube wall is ct = prk/hk, where p is the average ar- 
terial pressure. Since for large arteries p does not vary 
much, we take p constant. Then, we can study the be- 
haviour of the elastic energy stored in the walls, namely 
U = J{ay2E) dV = ipy2)J2^{Nkrllk/Ekhk). Using 
fk/hk constant, we find J7 ~ M in the case Ek = E. This 
proportionality between elastic energy and body mass 
can be understood as follows. At each heart beat, the 
energy imparted to the blood by the heart ||5| is given by 
(cardiac output)/ (heart rate) ~ M. Part of this energy 
is used to push the blood towards the capillaries, while 
another part is stored in the elastic arteries. The latter 
is needed to attenuate the pressure pulse, so that blood 
flow in the capillaries is smooth, allowing the exchange 
of nutrients and gases. Since this process is important 
irrespective of body size, hypothesis (vi) seems natural 
and is indeed observed in the case Ek = E. 

We now show that the same allometric exponent is 
found for MMR if we take the total elastic energy pro- 
portional to M instead of the mass of arteries. Assuming 
that the average arterial pressure has approximately the 
same value for all large arteries, we have a new auxiliary 
function, namely 

^ AE^'^h^'^ ^ 

F{rk,lk,Ekhk) = ^^-^ + XbY,Nkrllk + 

k=Q ^krk k=o 
N ^ w 3; 

Note that we may treat Ekhk as a single variable and that 
we have incorporated all constants in the Lagrange mul- 
tiphers. Requiring dF/dlj = dF/drj = dF/d{Ejhj) = 

gives us I3j — rij ^^^^ and Ejhj/rj — constant. In the 
capillary case, we again take the average pressure as con- 
stant, since the major pressure drop occurs in the arte- 
rioles. Using the impedance in Eq. (7) in the auxiliary 



this is equivalent to minimizing Z. Consider the to- 
tal impedance, Z — "^^^q Zk/Nk, for a network with 
no,ni, . . . ,nN-i fixed. For an organism of given mass 
M, we have to minimize with a given blood vol- 
ume Vf,, a given mass of thin- walled arteries (Ma ~ 
Nk^-KTklkhkPa, Pa IS the mass density of the arteries) 
and subject to the space-filling geometry. Defining the 
Lagrange multipliers (Ab,Afe,AQ), we minimize the auxil- 
iary function: 

^ AE^'^h^'^ ^ 
F{rk,lk,hk) = '^-J7l- + ^bY.^krllk + 

fe=0 NkTk k=0 

N N 

+ J2>'kNkll + XaJ2^krklkhk . (6) 

fc=0 k=0 

Here we have incorporated the constants in the Lagrange 
multipliers, and supposed that all vessels have the same 
mass density. To continue let us suppose that the Young's 
modulus is the same for all vessels {Ek = E). The re- 
quirements dF/dlj ~ dF/drj — dF/dhj ~ straigh- 

forwardly lead to (3j — n^ "^^^^ and hj/vj — constant. 
Using nj = n (hypothesis (v)), we obtain from Eq. (3) 
that the exponent for the maximum metabolic rate is 
b = 6/7 « 0.86. This result, larger than the basal 
value, is in very good agreement with experiment (see 
Table P). Additionally, we predict that capillary density 
scales as Nc/M ^ M"^/'' and maximum heart rate as 
B/M ~ A/^^/^, which is also in good agreement with 
data (see Table IJ). 

Until now we have considered only the large arteries. 
We must now consider the capillaries. Morgan and Kiely 
solved the Navier-Stokes equation for the liquid cou- 
pled to the Navier equation for the vessel wall in the limit 
of small radius. They found a highly attenued pressure 
wave with very small velocity c ~ {EkhkrkV/^/iujriy/^, 
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function we obtain (3j 



-1/3 



Therefore we have the 



-J - '-j 

same scenario as before and recover b = 6/7 « 0.86 for 
the MMR. 

The elastic properties of the vessels are important un- 
der heavy exercise. As observed by Womersley |21j, on 
the other hand, during basal metabolism the vessels can 
be approximated as almost rigid tubes because (a) Vb 
and Qo have their minimum functional values, and (b) 
heart rate is minimal, while the wavelength of the pres- 
sure pulse is maximal. (For example, in humans, where 
the aorta length = ~ 40 cm, Xbasai ~ 4 to and 
Xmax ~ 0, 9 to). The volume flow, defined in Eq. (4), 
can be written in terms of dp{z,t)/dz as 



Q ^ idp 
ujpa dz 



(9) 



If A >> /fe, we have that dp/dz « Apk/h, where Apk is 
the time-dependent pressure difference between the ex- 
tremities of the tube. Then the impedance is 



Poujlk 



(10) 



This is the same impedance found by Womersley 21] for 
a rigid tube with an oscillatory pressure gradient when 
is large and viscous effects are negligible. Now the local 
impedance coincides with the hydrodynamic resistance of 
the whole tube against Ap. To obtain /3, we must then 
minimize the following auxiliary function 



N 



N 



N 



k=0 



Nkri 



fe=0 



fe=0 



(11) 

The requirements dF/dvj — dF/dlj — furnish /3j ~ 
n ^^"^j characteristic of the so-called area-preserving 



regime. In the limit of small r^, viscosity dominantes 
so that the impedance is given by Poiseuille's formula 



(12) 



If we use this impedance in Eq. (11), we obtain Pj = 
rij . Therefore, assuming that Uj — n, Eq. (3) pro- 
vides 6 = 3/4 for basal metabolism. 

In summary, recent investigations have established 
that maximum metabolic rate scales with an exponent 
larger than that of basal metabolism. This suggests that 
basal and maximum metabolic rates are governed by dif- 
ferent mechanisms. We have shown, however, that both 
scaling regimes are predicted by a resource distribution 
network model 5] if we take into account the distinct 
hydrodynamic and elastic properties of the circulatory 
system at each metabolic state. The model correctly pre- 
dicts that MMR scales as M^/'', maximum heart rate as 
A/^^/'' and muscular capillary density as Af"^/''. The 
theoretical and experimental basis for the 3/4-law for 
basal metabolism were recently questioned by Dodds et 
al. However, Savage and coworkers extensively 

scrutinized the experimental data and confirmed that 
quarter-power scaling is pervasive in biology. Here we 
have shown that the 3/4-law is correctly predicted by 
the WBE model for mammals and birds when a differ- 
ent set of assumptions is adopted and the appropriate 
impedances are used. 
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